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Abstract 

Using Particle Image Velocimetry (PIV) in a new experimental plane Couette flow, we investigate 
the dynamics of turbulent patches invading formerly laminar flows. We evidence experimentally 
for the first time in this geometry the existence of large scale flows. These flows appear as soon as 
laminar and turbulent domains coexist. Spectral analysis is used to study the dynamical evolution 
of these large scales as well as that of the small scales associated with turbulence. We show that 
large-scale flows grow before turbulent spots develop and we point out the crucial role they play in 
the growth mechanism and possibly also in the emergence of organised patterns. 
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I. INTRODUCTION 


Transition to turbulence may occur through gradual supercritical scenarios as in Rayleigh- 
Benard convection or in Taylor-Couette flow when the outer cylinder is at rest. It is then 
well understood in terms of linear stability analysis. Beyond these cases, the transition to 
turbulence may also happen in a more sudden fashion that linear stability analysis fails to 
predict. Most wall-bounded shear flows fall into this last category and experience a subcrit- 
ical transition to turbulence involving the abrupt appearance of a disorganised state within 
the laminar phase. Plane Couette (PCF), Hagen-Poiseuille or plane Poiseuille flows are typ¬ 
ical examples of such a transition. The Reynolds number defined as Re = ULjp, where U 
and L are typical velocity and length scales and u is the fluid kinematic viscosity, being the 
natural control parameter of these systems, several critical values of Re can be pointed out. 
Following linear stability analysis, these flows should become turbulent beyond Re = Rec 
but experiments or numerical simulations show that turbulent patches can be sustained for 
values of Re exceeding some value Reg with Reg < Rec which is a global stability thresh¬ 
old below which any perturbation relaxes to the laminar state. In the PCF case, one can 
also define Ret beyond which the flow is always observed to be homogeneously turbulent 
and below which turbulence only fills a fraction of space Ft, hereafter called the turbulent 
fraction. It follows that a transitional range exists between Reg and Ret where laminar and 
turbulent domains coexist, where the flow is sensitive to finite size/amplitude perturbations 
and displays a strong hysteresis. 

In the nineties, two teams in Saclay and Stockholm have studied the transition to turbu¬ 
lence in PCF experimentally, pointing out its subcritical nature [H E] and characterising 
the growth of triggered turbulent patches thus denominated as spots [2113]. These turbu¬ 
lent spots consist of localised coherent structures made of small-scale longitudinal vortices 
and velocity streaks characteristic of wall turbulence and already evidenced and studied in 
boundary layer [1], channel flows |3] and numerical plane Couette flow in the pioneering work 
of Lundbladh and Johansson [H]. An other striking feature of PCF in the transitional range 
is the existence of a steady pattern of alternating turbulent and laminar stripes inclined 
with respect to the mean flow. These patterns were first observed experimentally by Coles 
and van Atta in Taylor-Couette flow[7l |8] and by Prigent et al. in PCF [9] and since then 
reproduced numerically by different teams in PCF[TDHI2] and channel flow [13]. 
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Spot growth has been carefully studied by the different teams cited above with similar re¬ 
sults: growth rates in the streamwise and spanwise directions have been found to be almost 
constant during the growth and to be increasing functions of the Reynolds number; waves 
at the spots’ edges have been observed to move at a speed slower than the spreading of 
turbulence. Regarding the growth mechanism, it could be due to stochastic nucleation of 
new streaks at the edge of the spot associated with a local destabilisation of the velocity 
profile at the laminar-turbulent interface, or to non-normal growth. A seemingly key element 
already pointed out in the work of Lundbladh and JohanssonfS] is the large-scale flow that 
develops around the turbulent growing spots with a kind of quadrupolar structure that has 
been observed numerically in PCF[6l[T3] and in model flows |15[[T6] but which is also present 
along the regular laminar/turbulent patterns [Hj. So far, they have never been observed ex¬ 
perimentally in PCF and only very recently in plane Poiseuille[I7j in spite of the important 
role they may play regarding spot growth and pattern sustainment as recently suggested 
by Duguet and Schlatter [Hj who proposed a mechanism combining stochastic nucleation of 
new streaks at the spot edges and their advection by large scale flows to explain the spot 
growth and its tendency to form oblique patterns simultaneously. 

In the present article, we aim at evidencing large-scale flows at the laminar-turbulent inter¬ 
faces in an experimental plane Couette flow and at quantifying them as much as possible all 
along the growth process. This would be a first step in the precise understanding of turbu¬ 
lence spreading in the transitional range of PCF. The article is organised as follows: we first 
describe our experimental set-up, the acquisition chain and the post-processing methods in 
im In section |III[ dedicated to quantitative measurements of large-scale flows in two very 
different laminar-turbulent coexistence situations, we particularly insist on the dynamical 
aspects. Implications of these results on spot growth mechanisms and also possibly on the 
onset of organised patterns are further discussed in [IV| 


II. EXPERIMENTAL SET-UP & METHODS 

A. Experimental set-up and measurement systems 

Plane Couette flow: the set-up sketched in figj^left was implemented to approach the 
ideal plane Couette flow developing between two infinite parallel plates. An endless 0.25 mm 
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thick plastic belt links two main cylinders. One of them is connected to a brush-less servo¬ 
motor which drives the system through a gear-reducer of ratio 9. This motor can achieve 
a maximal speed of 2500 rpm. Speed and acceleration are controlled with an accuracy 
of 2/1000. Four smaller cylinders guide the belt and enable to continuously adjust the gap 
between the two walls moving in opposite directions. This set-up is immersed in a tank filled 
with water. Water temperature is monitored by a thermo-couple with a 0.5° accuracy. The 
Reynolds number is defined as Re = Uh/u, where p is the kinematic viscosity of the fluid, U 
the belt velocity and h the half-gap between the two faces of the belt. Work presented here 
corresponds to a gap of 2h = 7.5 mm. Due to the experimental geometry, plane Couette 
profile is achieved within a. = 800 mm (along x) x = 400 mm (along z) area fitting 
our visualisation field of view. The corresponding aspect ratios are F^, = L^/h = 107 and 
Vz = Lz/h = 53. X, y and ^ and Uy and Uz are respectively the streamwise, wall-normal 
and spanwise directions and the associated velocities. The plane y = 0 corresponds to the 
half distance between the two plates. Unless otherwise specified, the tank is closed at the 
top by a transparent Plexiglas lid wherever the plane Couette profile is achieved in order 
to enable top and bottom symmetrical boundary conditions m- From various experiments 
not detailed here, critical Reynolds numbers have been obtained as Rg ~ 305 and Rt — 390, 
consistently with previous experiments and well-resolved direct numerical simulations. 

Velocity measurements: velocity measurements have been achieved by Particle Image 
Velocimetry (PIV). The Dantec system consists of a dual pulse Laser (Nd:YAG, 2 x 135 mJ, 
4 ns, 532 nm) and a CCD camera (FlowSenseEO, 4 Mpx). Image pairs are acquired at a rate 
of 5 Hz. The PIV Laser is mounted above the test section on a linear traverse which allows 
us to measure velocities over different y planes with a spatial resolution of about 1.3 mm 
along y. The laser sheet produced from a cylindrical lens gives us 2D velocity fields {Ux,Uz) 
in the streamwise/spanwise planes. The flow is seeded with particles coated with melamine 
and diameter {dp = 1 — 20 jim) that diffuses light at a shifted wavelength. A filter mounted 
on the camera lens is used to keep only wavelengths around this shifted wavelength, getting 
rid of unwanted light reflections on the vessel and plastic belt. The 1948 x 2048 pixels^ 
observation window corresponds to physical sizes of 196 x 208 mm^ or 305 x 310 mm^. An 
adaptive cross-correlation processing is applied to an initial interrogation area of 64 x 64 
pixels^ followed by a final interrogation area of 32 x 32 pixels^ with a 50% overlap. The 
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FIG. 1. Left: general sketch of the plane Couette flow set-up (a). Right: Ux space-time (t,z) 
diagrams obtained after a water jet has perturbed the flow (see text for details). U = 83.5 mm.s“^, 
Re — 330, 2h — 7.5 mm. Measurements are performed at y/h 0, y/h ~ 0.4 and y/h^ 0.5 from 
(b) to (d). Color online. 

corresponding spatial resolution is 2.45 mm, i.e. O.I6/1. Each instantaneous velocity field 
is further filtered by removing values exceeding 1.2?7. Resulting outliers are replaced using 
local median filters. 

Experimental protocol: our main goal is to focus on the coexistence of laminar and 
turbulent domains and on the associated large-scale flows. As explained in section |I| this 
situation occurs in the range of Reynolds numbers between Reg and Ret provided that a 
finite amplitude perturbation is added to the flow. Fig. [T]-right presents spatio-temporal 
diagrams of obtained after a water jet has perturbed the flow. In this case, the top lid is 
removed and water is blowed for 100 ms at 5 bars.The jet nozzle is situated at a: = 300 mm 
and .2 = 380 mm and is tilted by 45° with respect to the x direction. The flow is initially 
perturbed at large scales, but after a 10 s transient, the initial perturbation shrinks and 
a localized turbulent spot remains. From one realization to another, the fate of this spot 
drastically changes from decay to monotonic growth and even complex spatio-temporal 
dynamics. The three examples presented here illustrate the wide variety of dynamics that 
can be observed in the transitional range and call for more reproducible perturbations to 
study steady laminar-turbulent fronts. 
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To achieve such reproducible perturbations, we have used two very different protocols. The 
first relies on a permanent perturbation that consists of a bead positioned close to the 
y/h = —0.5 plane in the middle left of the PIV field of view (see bead position in fig. 
by means of a thin horizontal wire stretched between two vertical stalks placed between the 
guiding and the main cylinders. This type of permanent perturbation has already been used 
by Bottin and co-workers [18]. The wire diameter is (pwire/h = 0.1 and the bead diameter 
is (t^bead/h = 0.8. As explained in [12], a thin wire parallel to the streamwise direction does 
not influence the flow. As already mentioned by Bottin and co-workers, the bead diameter 
does not strongly affects the qualitative results and provided Re > Reg, no intermittency 
is observed but rather only growing spots are triggered. We studied the influence of the 
stalks on the transition to turbulence in our system when the bead is not present: their 
influence is hardly noticeable when considering critical Reynolds numbers and transition 
time to turbulence. When turbulence develops in a formerly laminar domain, it always 
occurs first around the bead. 

In an even simpler manner, a moving laminar-turbulent interface can be obtained by 
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FIG. 2. Snapshots obtained from our visualisation system. PIV field of view is within the red frame, 
(a) step experiment: turbulence invades the laminar domain from the left and right edges (along x 
direction) with a typical shape of an arrowhead, (b) bead experiment: a growing diamond patch of 
turbulence invades the laminar domain. The bead is situated at the orange dot. Visualisations are 
obtained from imaging an observation window of 800 x 400 mm^ fitting the area where the plane 
Couette profile is achieved. Flow is seeded with micron-size platelets (Iriodin) and illuminated in 
a y plane with a laser sheet produced by a rotating polygonal mirror. 

what we refer to hereafter as a step: a sudden increase of Re from a value below Rcg to a 
value above it. Practically, the belt is set into motion with an acceleration of 300 mm.s“^ 
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from Re = {) to Re = 450, i.e. above Ret so that the final state is expected to be fully 
turbulent. The same storyline always occurs: turbulence appears from the side edges as 
arrowhead shaped fronts as illustrated in fig. Actually, due to the experimental geometry 
(see fig. [^, the Reynolds number is higher between the guiding and the main cylinders so 
that turbulence is first triggered in these regions. For this type of experiments, the PIV 
field of view (red rectangle in fig. is shifted to the right of the set-up, hence we observe 
a turbulent front showing up from the right and slowly overrunning all the field of view by 
moving from right to left. 


B. Post-processing 


Measured quantities: to enlighten structures present in the flow, we use the autocorre¬ 
lation of Ux defined as: 


1 /‘+°° 

Ruu{x,Az,t) = — Ux{x,z + Az,t)Ux{x,z,t)dz (1) 

^0 J-oo 

with Rq = Ruu{x, 0, t). Ruu{x, Az, t) is further averaged over x and t to obtain i?„„(Az). 
Spatial spectral analysis is performed by computing Ux and L5, pre-multiplied power spectra 
obtained from Ux/z{kx,kz), the two-dimensional Ux and Uz Fourier transforms. The corre¬ 
sponding power spectra, Ex/z{kx,kz) = \Ux/zi.kx, kz)]"^ are azimuthally averaged to obtain 
Ex/z{k), with k = 2 t:/\ = \/k1 + k^. kx and kz being the wave numbers along x and .2 
respectively. Pre-multiplied power spectra are finally obtained as kEx/z{k). 

Validation on laminar and turbulent flows: in fig. Ux and Uz fields are shown in a 
(x,z) plane at Re = 450 (homogeneously turbulent case). These fields are the result of a 
sliding average of 5 snapshots acquired at 5Hz. Whatever the wall-normal plane considered, 
the velocity fields share the same attributes. An alternation of high and low speed streaks 
typical of turbulent wall shear flows are seen on Ux- Fluctuations are also observed on Uz, 
but their amplitude is twice lower and no obvious structure appears. plotted in fig. (c) 
exhibits a tip around 2h giving the characteristic size of the streaks as about four times the 
half-gap. This typical size is in line with the numerous previous estimations in wall-bounded 
shear flows. The time averaged pre-multiplied spectrum kEx drawn in figj^ (d) has a clear 
peak which corresponds to streaks energy. The associated wave number is coherent with the 
one found from the autocorrelation. The laminar flow spectrum is much lower and does not 
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(a) (b) (c) (d) 

FIG. 3. From left to right: velocity fields in the x — z plane, Ux (a) and Uz (b), U = 119,4 mm.s^, 
Re = 450, y/h 0. Antocorrelation of Ux along z averaged along the streamwise direction and in 
time (c) and ID time averaged power spectrum (d). Color online 

show any peak. for the same laminar flow (not shown) is that of a white noise giving 
the limitations of our PIV measurements. 


III. RESULTS 

A. Evidence of large scale flows 

For both perturbation protocols (steps and localized bead), the flow is fully laminar at 
the beginning of the experiment. Turbulence appears either at the streamwise edges of the 
experiment or around the bead as illustrated in fig. Then, it gradually invades more and 
more space before reaching a stationary turbulent fraction (see caption for details). This 
turbulent fraction may be equal to one provided the final Reynolds number is higher than 
Ret or less than one otherwise. In any case, the flow is most of the time in a state where 
turbulent and laminar areas coexist. Fig. [^presents premultiplied power spectra associated 
to Ux (a) and Uz (b) obtained from both types of experiments. The spectra are calculated 
at times corresponding to intermediate turbulent fractions and are averaged over 2 s (10 
consecutive samples). Two peaks are visible, one around A/h ~ 4 — 5 and the other around 
X/h ~ 40. The former corresponds to the streaks as discussed in section |n| the latter to 
large-scale structures. The peak associated with streaks is five times more intense on the 
streamwise component than on the spanwise one. This ratio is consistent with the case where 
the flow is homogeneously turbulent as discussed in|TIj Regarding the large-scale wavelength, 
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the same energy level is observed for both components. 

We have performed our measurements in various y/h planes. Regarding the streaks energy, 
the peak position and intensity are unaffected by the measurement plane while the energy 
distribution around the large wave number is. Large-scale peak associated to kE^ is more 
pronounced in planes closer to the belt. kE^ spectra are almost unaffected by the measure¬ 
ment position. To summarize, it is easier to detect large-scale flows on kE^ than on kE^., and 
for kEx, large-scale flows have a stronger signature on measurements performed on planes 
shifted from the midplane. 

The scale separation allows us to define without ambiguity a cut-off wavelength at Ac/h = 24 
to extract the large-scale and small-scale flows with 4*^ order Butterworth low-pass and high- 
pass filters. Fig. [^presents snapshots of and of the wall-normal vorticity (uy) associated 
with large and small scales. All snapshots correspond to the same given time along the 
growth of a spot around the bead at Re = 347 when the spot fills about half the PIV field 
of view. The low-pass filtered vorticity and velocity fields evidence a recirculation area all 
around the spot which corresponds to the quadrupolar flow obtained numerically by Duguet 
and Schlatter [Hj and Lagha and Manneville[T6]. The small-scale flow consists of long co¬ 
herent streamwise velocity streaks visible from the vorticity (f). These structures are much 
less coherent along x when considering Uz (d). These disorganized and spatially fluctuating 
structures are the footprint of turbulence inside the spot and thus allow to discriminate 
turbulent from quasi-laminar region. Considering fi^c and if together shows that the 
large-scale recirculation is present not only within the laminar area but also in the turbulent 
one. 

We have shown that along the spreading of turbulence into formerly laminar plane Couette 
flow large and small-scale flows coexist. Small-scale features are associated with streaks 
within the turbulent areas while the large-scale recirculation spreads over laminar and tur¬ 
bulent areas. This recirculation is linked to the laminar-turbulent coexistence and is involved 


in the spreading of turbulence as discussed in section |lVj The temporal evolution of these 
structures is examined in the following sections. 
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B. Growth dynamics 


Fig. shows [4 fields at different times along a step experiment during which a turbulent 
front invades the laminar flow from right to left in the pictures. In the first snapshot, the 
flow is laminar and the [4 field is almost uniform. In the second one, a narrow turbulent 
area is already visible at the top right corner and a large scale recirculation very similar to 
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FIG. 4. Top: dimensionless 14 (a) and Uz (b) premultiplied power spectra for typical step (trian¬ 
gles) and bead (circle) experiments. Measurement are done dX yjh 0.2 and Re = 450 for the step 
and at ^//i ^ 0 and Re = 403 for the bead. Bottom: velocity and vorticity field in (x, z) plane for 
a bead experiment. Uz at large (c) and small (d) scale, lOuy at large (e) and ujy at small (f) scale. 
Measurements are done dX yjh 0, Re = 347. The bead is situated at the black dot. A dashed line 
corresponding to the turbulent spot contour obtained from (f) is shown for all snapshots. Color 
online. 
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that observed around the turbulent spot in fig. |4]-c spreads over half the picture. In the third 
snapshot, this large-scale structure has moved to the left along with the laminar-turbulent 
front. On the last snapshot the recirculation is present only close to the laminar-turbulent 
interface, while deep inside the turbulent area, far enough from the front, the velocity field is 
as in fig. 1^ where the flow is homogeneously turbulent. This exemplifies the fact that large- 
scale flows go along with the growth dynamics of turbulent areas. To better understand this 
interplay, we now consider the time evolution of the premultiplied power spectra. 

Fig. 1^ show spectrograms associated with a step experiment (a) and a bead experiment (b). 
In both cases, the scenario is similar: during a first phase (1100 h/U for the step, 300 h/U 
for the bead), spectra are almost uniform, do not present any peak and their magnitude are 
low: the flow is laminar. The noise visible at large scales defines the resolution of our PIV 
setup. During a second phase, energy emerges at large wave numbers (greater than lOh), 
quickly followed by a third phase where energy raises around A/h ~ 4 — 5 while the energy 
at large scales still grows. This corresponds to the emergence of large-scale flows and to 
the establishment of well organised streamwise streaks {i.e. turbulence) on larger and larger 
areas. Note that the corresponding energy level are two to three times greater than the 
noise present in the first phase. In a final phase corresponding to a state where turbulence 
has invaded the whole flow, the large scales are damped and all the energy gathers around 
A/h ~ 4 — 5. We insist on the fact that the energy associated with large-scale structures 
appears earlier than the energy finked to turbulence. Spectra associated to five instants 
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FIG. 5. Uz snapshot during a step experiment, (a) tU/h — Q2 , (b) tU/h — 370, (c) tU/h — 551, 
(d) tU/h — 933), Re=450. Vectors correspond to the {Ux,Uz) large scale velocity field. Color 
online. 
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FIG. 6. Left: spectrograms of the streamwise pre-multiplied energy kE^. measured at y/h ^ 0.2. 


Right: spectra extracted from the left figures at instants marked by grey fines, a-c: step at 


y/h 0.2 and Re = 450 , b-d: bead at y/h ^ 0 and Re = 403. Color online 


representing these three phases are extracted and plotted on fig. |^c-d (see caption for 
details). On these plots, one better appreciates the low level of the noise during the laminar 
phase and the hierarchy of energy level of large and small-scales as time passes. Large scales 
are associated with the peak growing at X/h = 40. The amplitude of both peaks are similar 
even if the large-scale one is somehow higher. 

We now focus only on the spectral components associated to large and small-scales. The 
small scale {i.e. the streaks) component Astr is defined as the average of kE^ over 4 < 
X/h < 5 and the large-scale component Alsf as the value of kEz at X/h = 40. Their time 
evolutions are presented on fig. [^a and fig. 0b respectively. Each curve corresponds to a 
phase averaging of five realisations at a given value of Re. Before performing the phase 
average, a time sliding average over 1 second (5 samples) is applied. Regarding streaks, the 
amplitude Astr remains roughly constant at a very low level before suddenly increasing 
during about 800 h/U to reach a well defined plateau. The level of this plateau (shown as 
horizontal lines in the figure) increases with the Reynolds number. The amplitude Alsf 
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time evolution differs more from one value of Re to another but three tendencies still appear, 
(i) For Re ^ 330 Ai^sf increases monotonously to reach a plateau at relatively low level: 
large-scale flows are present but the growth is limited and the turbulent spot will not invade 
the whole spanwise extent of the experiment. These spots survive rather than grow. By 
contrast, if Re > 330, A^sf increases rapidly to reach a maximum. Then it either remains 
around this maximum (ii) or relaxes until stabilising on a more or less well defined plateau 
(hi). Distinction between (ii) and (iii) corresponds to spots eventually hitting the spanwise 
boundaries {Re > 350) or not. Case (iii) includes situations with Re > Ret for which 
Alsf vanishes rather than relaxes. Note that the phase averaging may hide an important 
dispersion between various realisations obtained at the same Re. Cases (ii) and (iii) are 
better understood in fig. [^c and 0cl where large and small-scale amplitudes are plotted 
together for two given realizations at Re = 340 and Re = 403 respectively. For Re = 340, 
from t = 0, large-scale amplitude starts increasing before streaks amplitude does. Then 
both grow together before reaching plateaus. For Re = 403, at a time around 10 h/U, 
both amplitudes start to increase but this increase is faster for the large-scale than for the 
streaks. A^sf stabilizes to its maximal value while the streaks one keeps increasing. Once 
Astr reaches its plateau, A^sf drops sharply to eventually vanish. This corresponds to 
situations above Ret for which the ultimate state is homogeneously turbulent and large- 
scale flows are thus transient. Note that an overshoot is present at large scale. For range of 
Re below Ret but above Re ~ 350) the spot develops until it reaches the spanwise edges of 
the experiment. Then, the turbulent fraction remains almost steady but the turbulent spot 
reorganises as oblique bands. Large-scale flows are present all along this reorganisation (the 
plateau mentioned above) but also once the bands are stabilised as shown numerically by 
Duguet and Schlatter [Ill- 

Fig. 1^ shows the evolution with the Reynolds number of the maximal amplitude associated 
to the large scales (defined as the maximum of the curves presented in fig. [^a and of the 
streaks amplitude (defined as the value of the plateau in fig. 0b- Error bars displayed 
in this figure correspond to the dispersion between various experiments performed at the 
same Re we alluded to above. If the streaks amplitude is increasing smoothly in the range 
of Re considered, a cross-over between two regimes may nevertheless be identified around 
Re ~ 350. At these values of Re, the behaviour of the large-scale maximal amplitudes 
changes abruptly from steep to moderate increase. All this would discriminate between the 
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case (ii) and (iii) described above. This point will be further discussed in the next section. 


IV. DISCUSSION AND CONCLUSIONS 

PIV measurements have been performed in plane Couette flow around laminar-turbulent 
interfaces associated to growing turbulent patches. Two very different protocols have been 
used, one with a localised perturbation achieved by a bead placed in the center of the 
channel and one with no specific perturbation but with turbulence generated at entrance 
denominated as step. In both cases, turbulence grows provided Re is higher than Rcg. 
In the bead case, turbulence consist of a localised diamond like spot growing around the 




(a) (b) 




(c) (d) 

FIG. 7. Time evolution of streaks Astr (a) and large-scale Alsf (b) amplitudes for several 
Reynolds numbers. Time evolution of streaks Astr{^^A) and large-scale flows ^^.^'^^^(blue) ampli¬ 
tudes for given realizations at Re = 340 (c) and Re = 403 (d), y/h ^ 0. 
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FIG. 8. Evolution with the Reynolds number of the maximal amplitude associated with the large 
scales (blue) and of the streaks plateau amplitude (red), y/h ~ 0. 

bead while in the step case turbulence arises from the stream-wise edges as large arrow 
shaped fronts. Spatial spectral analysis has been used to evidence for the first time the 
existence of large scale recirculations at the laminar-turbulent interface in an experimental 
plane Couette flow. Existence of large-scale flows in situations where laminar and turbulent 
domains coexist has already been explained by flow rate conservation arguments in presence 
of the overhang existing at the laminar-turbulent boundary [g Ej or through additional 
pressure fields induced by the distribution of Reynolds stresses [IS] . 

As soon as laminar and turbulent areas coexist in the experiment, even if turbulence 
eventually decays, large scale flows develops not only around the fronts but also far into 
the laminar and turbulent domains. Around the localised spots, they consist of quadrupolar 
recirculation cells as already reported in numerical simulations of plane Couette flow miiisi 
and idealised shear flow na but also experimentally in plane Poiseuille na. Quantitative 
comparison to previous numerical works is difficult since the authors did not provide actual 
values to estimate the large-scale flow intensities. Nevertheless, we have performed numer¬ 
ically a similar but less exhaustive study on ChannelFlow [22l [23] where we initiate spot 
growth with the initial perturbation used in [6] . The results obtained so far are in good quan¬ 
titative and qualitative agreement with the present work, whenever the intrinsic differences 
between experimental and numerical setup do not prevent such comparison. Regarding our 
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work, large-scale flow intensities are similar between the two types of experiments we have 
considered but they strongly depend on the Reynolds number as evidenced in flg. and 
on the turbulent fraction. This dependence on the turbulent fraction can be understood in 
the framework of flow rate conservation argument evoked above: the larger the turbulent 
fraction, the longer the overhang and thus the larger the flow rate imbalance that has to be 
compensated by stronger large scale flows. 

According to our understanding and experimental observations, large-scale flows exist all 
along the growth process till turbulence has spread over the whole spanwise direction. Then 
two scenarios are possible: if Re < 330, large-scale flow exist but spots survive rather than 
grow. If spots ultimately grow {i.e. if Re > 330) large-scale flows show rapid increase 
and eventually stabilise on finite value associated to a steady turbulent fraction. Our re¬ 
sults suggest an additional distinction between experiments performed at Re below or above 
Re ~ 350. In the latter case, triggered spots reach the spanwise extend of the experiment 
before reorganising as steady regular patterns. In the former, the steady ultimate turbulent 
fraction is too small to accommodate such patterns and the triggered spot starts reorganising 
as inclined patterns before reaching the spanwise edges. The final state is a wandering mix 
of inclined fractioned stripes whose turbulent fraction highly fluctuates around a stationary 
value. This “critical” Re ~ 350 is consistent with Prigent’s work who reported that organ¬ 
ised steady patterns can be observed only above Re ~ 340 while below stripes are not stable 
and one mostly observes stretched spots|20]- Some preliminary numerical simulations similar 
to the present study confirm most of these observations and points out the resurgence of the 
large scales when the turbulent spot later reorganises to form laminar/turbulent stripes as 
those observed in large domains [MI] .An other intriguing and interesting feature is the role 
played by large-scale recirculations in the spreading of turbulence. On the one hand, we have 
shown that these recirculations are present as soon as laminar and turbulent domain coexist, 
growth being observed or not, but on the other hand, Schumacher and Eckhardt [I3| attest 
that no growth is observed when what they call the outflow is suppressed around a localised 
turbulent spot in a shear flow very similar to plane Couette. In a streamwisely constrained 
plane Couette geometry where no large scale flow is present, Duguet et al. ini evidence a 
purely stochastic growth of turbulent spots but at rates far below what is observed in more 
realistic situations. Large-scale flows thus appear as a key ingredient of the growth even 
if they are not sufficient. Another point of debate is whether the large-scales are induced 
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by or induce the growth. As seen in fig. large-scale amplitude increases before streaks 
amplitude does, implying that growth is initiated by these large-scale fiows. How do they 
actually act on the growth, by simply advecting the turbulent region or by locally modifying 
the laminar profile to initiate a growth by destabilisation jl], remains an open question. The 
energy transfers between large and small-scales evidenced in the present article should be 
scrutinised. We are currently studying the growth rates of turbulent spots in bead exper¬ 
iments through high resolution visualisations in an attempt to refine our knowledge of the 
spot growth mechanisms. We hope to be able to discriminate the role played by large-scale 
flows in this process. 


V. ACKNOWLEDGEMENTS 

The authors would like to thank Y. Duguet, P. Manneville and O. Cadot for helpful 
discussions. The KTH and Saclay groups are acknowledged for sharing their experience 
regarding the experimental set-up, particularly F. Daviaud for lending the polygonal mirror. 
The master students L. Colo, G. Anyigba, E. Bach, C. Jacob, E. Sieffert, P. Yang and 
Y. Zhu are also acknowledged for their contributions. E. Moisy is acknowledged for is 
PIVMat toolbox. This project was supported by ANR Jeunes chercheurs/Jeunes chercheuses 
QANCOUET. 


[1] F. Daviaud, J. Hegseth, and P. Berge. Subcritical transition to turbulence in plane couette 
flow. Physical Review Letters, 69:2511, 1992. 

[2] N. Tillmark and P. H. Alfredsson. Experiments on transition in plane Couette flow. Journal 
of Fluid Mechanics, 235:89, 1992. 

[3] O. Dauchot and F. Daviaud. Finite amplitude perturbation and spots growth mechanism in 
plane Couette flow. Physics of Fluids, 7:335, 1995. 

[4] M. Gad-El-Hak, R. F. Blackwelderf, and J. Riley. On the growth of turbulent regions in 
laminar boundary layers. Journal of Fluid Mechanics, 110:73, 1981. 

[5] D. R. Carlson, S. E. Widnall, and M. F. Peelers. A flow-visualization study of transition in 
plane Poiseuille flow. Journal of Fluid Mechanics, 121:487, 1982. 


17 



[6] A. Lundbladh and A. V. Johansson. Direct simulation of turbulent spots in plane couette 
flow. Journal of Fluid Mechanics^ 229:499, 1991. 

[7] D. Coles. Transition in circular Couette flow. Journal of Fluid Mechanics^ 21:385, 1965. 

[8] CW Van Atta. Exploratory measurements in spiral turbulence. Journal of Fluid Mechanics^ 
25:495, 1966. 

[9] A. Prigent, G. Gregoire, H. Chate, O. Dauchot, and W. Van Saarloos. Large-Scale Finite- 
Wavelength Modulation within Turbulent Shear Flows. Physical Review Letters^ 89:014501, 
2002 . 

[10] D. Barkley and L. S. Tuckerman. Computational Study of Turbulent Laminar Patterns in 
Couette Flow. Physical Review Letters^ 94:014502, 2005. 

[11] Y. Duguet, P. Schlatter, and D. S. Henningson. Formation of turbulent patterns near the 
onset of transition in plane Couette flow. Journal of Fluid Mechanics^ 650:119, 2010. 

[12] J. Rolland and P. Manneville. Ginzburg-landau description of laminar-turbulent oblique band 
formation in transitional plane couette flow. The European Physical Journal B-Condensed 
Matter and Complex Systems^ 80:529, 2011. 

[13] T. Tsukahara, Y. Kawaguchi, H. Kawamura, N. Tillmark, and P. H. Alfredsson. Turbu¬ 
lence stripe in transitional channel flow with/without system rotation. In Seventh lUTAM 
Symposium on Laminar-Turbulent Transition^ page 421, 2010. 

[14] Y. Duguet and P. Schlatter. Oblique laminar-turbulent interfaces in plane shear flows. Physical 
Review Letters^ 110:034502, 2013. 

[15] J. Schumacher and B. Eckhardt. Evolution of turbulent spots in a parallel shear flow. Physical 
Review F, 63:046307, 2001. 

[16] M. Lagha and P. Manneville. Modeling of plane Couette flow. 1. Large scale flow around 
turbulent spots. Physics of Fluids^ 19:094105, 2007. 

[17] G. Lemoult, J.-L. Aider, and J.-E. Wesfreid. Turbulent spots in a channel: large-scale flow 
and self-sustainability. Journal of Fluid Mechanics^ 731:R1, 2013. 

[18] S. Bottin, O. Dauchot, and F. Daviaud. Intermittency in a locally forced plane couette flow. 
Physical Review Letters^ 79(22):4377, 1997. 

[19] S. Bottin, O. Dauchot, F. Daviaud, and P. Manneville. Experimental evidence of streamwise 
vortices as finite amplitude solutions in transitional plane Couette flow. Physics of Fluids^ 
10:2597, 1998. 


18 



[20] A. Prigent. La spirale turbulente: motif de grande longueur d’onde dans les eeoulements 
eisailles turbulents. PhD thesis, Universite Paris Sud-Paris XI, 2001. 

[21] Y. Duguet, O. Le Maitre, and P. Schlatter. Stochastic and deterministic motion of a laminar- 
turbulent front in a spanwisely extended Couette flow. Physieal Review 84:066315, 2011. 

[22] J.F. Gibson. ”Channelflow: A spectral Navier-Stokes simulator in C++”. J. F. Gibson, U. 
New Hampshire, 2014, channelflow.org. 

[23] J.F. Gibson, J. Halcrow, P. Cvitanovic. Visualizing the geometry of state space in plane 
Couette flow. Journal of Fluid Meehanies, 611:107, 2008. 

[24] Before having installed the top Plexiglas lid, the boundary condition at the top was a free 
surface and we observed slow global toward the free surface. Since the experiment has been 
closed at both top and bottom, we do not observe any asymmetry in any of our results. 


19 



